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Abstract 

We present an exact formula for moments and cumulants of several 
real compound Wishart matrices in terms of an Euler characteristic 
expansion, similar to the genus expansion for complex random matri- 
ces. We consider their asymptotic values in the large matrix limit: 
as in a genus expansion, the terms which survive in the large matrix 
limit are those with the greatest Euler characteristic, that is, either 
spheres or collections of spheres. This topological construction moti- 
vates an algebraic expression for the moments and cumulants in terms 
of the symmetric group. We examine the combinatorial properties dis- 
tinguishing the leading order terms. By considering higher cumulants, 
we give a central limit-type theorem for the asymptotic distribution 
around the expected value. 

1 Introduction 

Wishart matrices are random matrices of the form AX T BXC, where X is a 
matrix with Gaussian entries and A, B and C are constant matrices which 
may be required to be symmetric, positive definite, or the identity matrix. 
This matrix ensemble was introduced in Wishart's paper of 1928 [18] in 
order to study the expected covariances of multivariate data. The moments 
of Wishart matrices, as well as more general symmetric polynomials in their 
eigenvalues, are studied in such papers as [3,5,6,9, 10]. The behaviour of 
these quantities in large matrices and the connections to free probability are 
considered in papers such as [5,6,15]. 

In this paper, we exhibit an exact formula for the expected value of 
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E (TV (X T D 1 XD 2 • • • X T D 2ni -iXD 2ni ) ■ ■ ■ 

TV (X T D2n-2n r +lX D2n-2n r +2 ■ ■ ■ X T D2n~\XD2n) ) 

= ^TV 7 - 1 - 7+/2 (D 1 ,..., J D 2n ) 

TV 

where the sum on the right hand side is over all pairings tt on 2n elements, 
interpreted as permutations, ix is a permutation derived from tt, and 7_ and 
7 + are fixed permutations dependent on the left hand side expression. If we 
wish to calculate cumulants, rather than moments, we sum only over n such 
that the subgroup (74.7Z ,7r) is transitive. 

If the permutations are interpreted as the information to construct a 
topological cell-complex, this expression becomes a genus expansion similar 
to those for complex matrix ensembles (see [14, 15, 19]), but which also in- 
cludes terms corresponding to non-orientable surfaces. As in the complex 
cases, the order of the terms depends only on the Euler characteristic of 
the corresponding surface, allowing us to distinguish which terms survive 
in the asymptotic limit. The leading terms are again those with the high- 
est Euler characteristic. The terms contributing to the mixed moments of 
traces are arbitrary surfaces, so the highest order terms are those whose 
corresponding surface consists of as many disjoint spheres as possible. The 
terms contributing to cumulants of the traces are connected surfaces, so the 
leading terms are those corresponding to spheres. 

This topological interpretation is also conducive to studying several in- 
dependent Wishart matrices, or the more general case of Wishart matrices 
indexed by a Hilbert space whose covariances are given by the inner product. 
This construction also allows us to handle expressions in which the random 
matrix X and its transpose appear arbitrarily in the products of matrices, 
rather than alternatingly, as in the usual Wishart matrix model. 

The topological interpretation motivates the algebraic expression for 
these expected values, following the expression in [15] for the complex case 
based on the interpretation of a map or hypermap as a pair of permutations 
in [7]. We now must consider non-orientable surfaces as well as orientable 
ones. The algebraic formulas in [7] require a consistent orientation on the 
surface, so we construct orientable two-sheeted covering spaces of our origi- 
nal surfaces and extend the algebraic formula to this situation. 

In Section [2] we will outline the notation and results we will be using and 
define the matrix model which we will be considering. To simplify our no- 
tation, we will delay discussion of the several matrix model and expressions 
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with arbitrary transposes until Section [7j however, most of the arguments 
are identical. In Section [3] we provide the geometric motivation for the alge- 
braic formulas using a construction similar to [14, 19]. For a rigorous proof, 
we prefer to use an algebraic construction for the moments and cumulants 
using the representations for graphs in surfaces developed in [7] and used in 
the complex Wishart case in [15], which we describe in Section |4j The proof 
of the main formula is given in Section [5] and its asymptotic behaviour is 
discussed in Section [6] 

2 Notation, Definitions and Lemmas 

Notation 2.1. Denote the set {l,...,n} by [n] and the set {l,...,n} U 
{-l,...,-n} by [±n). 

We write the usual trace Tr and the normalised trace tr := j?Tr. 
2.1 The Matrix Model 

Definition 2.2. We will denote by X an M x N matrix whose ijth entry is 
^j^fijj where the are independent real iV(0, 1) random variables (that 
is, centred Gaussians with standard deviation 1). 

A matrix of the form X T D^XDf,^, where is an arbitrary constant 
M x M matrix and Dk+i is an arbitrary N x N matrix, is called a real 
Wishart matrix. 

We will be considering the statistics of these matrices as N —* oo, and 
we assume M — > oo at the same time. We will continue to consider the 
normalised trace Tr, even when we are considering an M x M matrix. We 
will also assume that the large N limits of the normalised traces of products 
of Dk matrices approach a finite limit. 

Let ni, . . . ,n r be positive integers. We will be considering the moments 
and cumulants of the traces of products of Wishart matrices, that is random 
variables of the form 

(T 
X D 2(n 1 +-+n k _ 1 )+lXD 2 ( ni+ ... +nk _ 1 ) + 2 ' ' ' 

X T D 2{n 1 +-+n k )-l X D 2 ( ni+ ... +nk )) ■ (1) 

There are no requirements on the matrices: they need not be symmet- 
ric, and none are required to be equal. Both the subscripts of the matrices 

-Dfc and the appearances of transposes appear naturally in our formula. We 

(ft) 

will denote the transpose of by _D_^ and the ijth entry of by oL . 
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Since we only consider the traces of products of matrices, expressions 
which begin with a constant matrix can be calculated by cycling any initial 
matrix to the end of the expression inside the trace. 

2.2 Partitions 

Definition 2.3. A partition of a set S is a set of subsets V±, . . . , V rn C S 
called blocks such that V* 7^ for 1 < % < m, Vi fl Vj = for alH 7^ j, and 
V\ U ■ ■ • U V m = S. We denote the set of all partitions of S by V (S), and 
the set of all partitions of [n] by V (n) . 

A partitions whose blocks contain exactly two elements is called a pair- 
ing. We denote the set of all pairings of a set S by V2 (S), and the set of all 
pairings of [n] by V2 (n)- 

2.3 Moments and Cumulants 

Definition 2.4. The nth mixed moment of random variables X±, . . . , X n is 
defined as 

a n (X 1 ,...,X n ) :=E(AV--X„). 

Definition 2.5. Cumulants are functions k±, k2, ■ ■ ■ defined so that k n is an 
n-linear function, and satisfying the moment-cumulant formula: 

a n (X 1 ,...,X n )= II k m (X h ,...,X im ) (2) 

vreP(n) V={h,...,im}en 

for all random variables X\, . . . , X n and all n > 0. 

It can easily be shown by induction that for all n, the moment cumulant 
formulas up to n along with the cumulants k±, . . . , fc n _i uniquely define the 
multilinear function k n . The first cumulant is the expected value, and the 
second is the covariance. 

We denote the term on the left hand side of the moment-cumulant for- 
mula corresponding to a given partition tt G V (n) by 

k w (Xi, . . . , X n ) := J J k m (Xjj , . . . , Xi m ) . 

V={h,...,i m }en 

2.4 The Wick Formula 

The following formula allows us to express the expectations of products of 
Gaussian random variables in terms of a combinatorial construction. 
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Lemma 2.6 (The Wick Formula). Let {f\ : A G A} be a centred family of 
Gaussian random variables (that is, each f\ is a linear combination of ran- 
dom variables from a collection of independent N (0,1) random variables). 
Then 

E(/a 1 ---/aJ= £ II E(/ A J Ai ). 
Tev 2 (n) {k,i}ew 

In particular, if n is odd, there are no pairings, and the expected value 
is 0. 

For a proof, see [13], theorem 1.23. 



2.5 Traces Along Permutations 

For a permutation it G S n with cycle notation (c±, . . . , c ni ) ■ ■ ■ (c ni +...+ rv _ 1 +i, 
. . . , c„), we define the trace along ir of matrices Ai, . . . , A n as 

1% (A,,..., A n ) := Tr(A Cl • • • A Cni ) • • • Tr(^ Cni+ ... +rir _ i+1 • • • A Cn ). 

Because the trace is cyclic, this expression is well defined. For permutations 
on a set of signed integers, we let a negative integer denote a transpose: if we 
let A { p = A k while A { ~^ = A\, and if ir = (c 1 , . . . ,c n J • • • (c ni+ ... +rir _ 1+1 ■ ■ ■ 
c n ) is a permutation of a subset of [±m], then 

Tr.^,...,^) :=T1^»...aW^ 

Tr ( ^( sgn ( Cll i+-+™--i +1 )) . . . ^(sgn(c n )) 

V | c n 1 +---+n r _ 1 +l| M 

We will be making use of the following folklore result: 
Lemma 2.7. Let Ai,...,A n be matrices such that A k is an x N a ^y.) 

(k) 

matrix. If we denote the ijth entry of matrix A k by a- , then 
Tr CT (A i,..., A n )= Y of] ■■■oS n ] . 

CT V ' U > «V(1) *nV(n) 

i<n<ATi 

l<i„'<JV n 

This formula may be proven by standard calculation. The subscripts on 
the matrices Ak may be replaced by the elements of a finite index set and a 
by a permutation of that set, allowing us to use this formula in the situation 
where the index set is a subset of [±ra] . 
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3 Geometric Motivation 



In this section we sketch a method of calculation using surface gluings. We 
delay a rigorous proof until Section [5] which will use only Lemma 3.2 from 
this section. 

We begin by examining the moments of the random variables defined 
in Q. We will let n = ri\ + • • • + n r {n^ as in Q). 

Notation 3.1. Let 7 := (1, . . . , 2n\) • • • (2 {n\ H h n r _i) + 1, . . . , 2n) £ 

Lemma 3.2. ^4 mixed moment of the Y^ may be expressed in terms of the 
entries of the matrices: 

a n (Yi, ...,Y r ) 

= N~ n - r V V d {1) d (2) ■ ■ ■ d {2n ~ 1] d {2n) 

~" L / '1' 7 (1) J2i 7 (2) i 2n~l* 7 (2n-l) hnjy(2n) 

ireV2 l<il,...,i2n<N 

i<h,-,h n <M 

II MfaikM- (3) 
{fc,i}e7r 

Furthermore, since we can express the product of expected values as: 

n w(f , \ f I3 i k = k and j k = ji for all {k,l} £ n 
^UikikJim) \ o otherwise { > 



we may replace the expected value expression with a condition on the indices: 
a r (Y 1 ,...,Y r ) = N~ n - r V V df] 41 ••• 

7re"P2(2n) i:[2n]->[Af]:i=io7r 
j:[2n]-»[M]:j=jfo7T 

d (2n-l) d (2n) 

«2n-li 7 (2n-l) 32njj(2n) 

where i (k) = ik and j (k) = jk- 
Proof. We can write the moment: 

a r (Y lr -- ,Y r ) = E(yi-..y r ) 

= E (tr (X T D 1 XD 2 ■ ■ ■ X T D 2ni ^XD 2ni ) ■ ■ ■ 

(T 
X D 2(n 1 +-+n r -i)+l XD 2(n 1 +-n r - 1 )+2 ' ' ' 

X T D 2n ^XD 2 n)). 
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We can interpret this expression as a trace along permutation 7, so according 



to Lemma 2.7 



■(«,•■■, Y r ) = E (N-'Ti-, (X T D U XD 2 , X T D 2 „- U X, D 2 „)) 



N- V E(xJ lh ^l m X hj ^ m 

l<i 1 ,...,i 2n <N 

l<31,-,32n<M 

-11 ) 

N- n ~ r V d (2) ••• 

l<jl,...,i 2 „<A r 
l<il,-j2n<M 

j(2»-l) ,(2n) E ,, , v 

%n-lS( 2 n-l) U j2ni 7 (2n) ^ WHJl J I2n32n ) " 



Applying the Wick formula (Lemma 2.6) to the expected value expression 
in the above, we get: 

a r (y l5 ...,y r ) = N- n - r V 4 2 !- ••• 

' y 11 ' 7 > Z-^ '1*7(1) 3l3 1 {2) 

l<il,...,i 2 n<iV 

i<ii,-,i2u<M 

>(2n-l) ,(2n) y> TT E , f . . f . . n 

n£V 2 {k,l}eir 

= N- n - r V V d {2) ■ ■ ■ 

iy %H(i) a j2j 7(2) 

7reP 2 l<u,...,j 2n <iV 
l<ilvj2n<M 

j(2n-l) ,(2n) TT TO f /•. . f..) 

%n-lS(2n-l) U J2nJ 7 (2„) 11 ^ \J WkJ H3l ) 
{k,l}GlT 

where we have reversed the order of summation to express the moment as 
a sum over pairings. This proves ([3]). 

In the term corresponding to a given pairing it, the factor E (fi k j k fi l j l ), 
and therefore the entire term, will be equal to zero if fi k j k is independent 
from fi l j l . So we find: 

n , s _ / 1, %k = k and = ji for all {fc, /} G 7r 

^Vidkhm) "I o, otherwise 

proving @. We can interpret the conditions for a nonzero expected value 
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Figure 1: The faces constructed in the calculation of 
E (tr (X T D 1 XD 2 ■ ■ ■ X T D 5 XD 6 ) tr (X T D 7 XD 8 X T D 9 XD 10 )) , with 
the random variables fi 2 j 2 and fi 7 j 7 (from X and X T respectively) paired 
with each other, and fi 1 j 1 and fi g j g (both from X T ) paired with each other. 
The resulting edge identifications are shown. 

as a further set of constraints on the indices: 
a r (Y 1 ,...,Y r ) = N- n - r V V d 1 } df] ■■■ 

TreV 2 (2n) i:[2n]-t[N]:i=iam 
j:[2n]-*[M]:j=jon 

d (2n-l) d (2n) 

*2n-l« 7 (2n-l) hnj 1 (2n)' 

This completes the proof of the lemma. □ 

These constraints can be expressed geometrically as follows. For each 
random variable Yp., we construct a 2nfc-gon. We associate each edge with a 
random matrix term and each vertex with a constant matrix term, cyclically 
(counterclockwise) in the order they appear in the expression for Y^. We 
label each end of each edge with the index the edge and the vertex at that 
end share. An example is shown in Figure [TJ 

A pairing ir constrains i/% = i\ and jk = ji for each {k, 1} £ tt. We 
represent this by identifying the edges X^ and X{ so that the edge-ends 
which have been labelled if. and i\ line up, as do those labelled j\. and 
ji. Figure [T] shows two examples. A pairing of a random variable from 
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Figure 2: A vertex that appears in the surface resulting from the pairing 
{{1,9} {2, 7} {3, 4} {5, 10} {6, 8}}. To perform the twisted edge identifica- 
tions, the corners of faces containing the matrices D§ and Dg were nipped 
over, so the matrices D§ and Dg appear instead, since it is D§ whose 
iei^th entry appears, and likewise with Dg. The indices are identified as 
they would be in the expression for the trace Tr (DiD^D^Dg) . 

an X term with one from an X T term will result in an untwisted edge 
identification while a pairing of random variables from two X terms or two 
X T terms will result in a twisted pairing. The resulting surface thus may 
not be orientable. However, since each edge is paired with exactly one other, 
it will be a compact surface without boundary. 

To determine the constraints on the indices of the constant matrices, 
we examine the vertices of the resulting surface. Figure [2] shows a an ex- 
ample. Due to the presence of twisted edge identifications, some corners 
may appear flipped over relative to others. We replace their matrices with 
their transposes so the indices of the matrix as written always appear in 
counterclockwise order. 

We note that around a vertex the indices of the matrices are identified as 
they would be in the expression for a trace. When we express the moment 
in terms of normalised traces, each trace gives us a factor of N. The order 
of the term corresponding to a given pairing whose surface has j^V vertices 
will be N& v ~ n ~ r , or, since the number of edges (after identifications) is n 
and the number of faces is r, m~ 2r , where x = #V-#E + #F is the Euler 
characteristic of the surface. The highest order terms will be those with the 
greatest Euler characteristic. 

Remark 3.3. We note that this construction appears implicitly in many of 
the apparently different approaches to the calculation of moments. 

The face gluings we have constructed can be interpreted as the duals of 



9 



the Mobius graphs in [4,16]. 

In the constructions used in [3,9, 10], two graphs are constructed whose 
union consists of cycles whose lengths are 2m, . . . , 2n r . We can interpret 
these cycles as the faces in our construction, although the vertices are edges 
and the edges are vertices. The moment is then expressed as a sum over all 
pairings of the vertices, equivalent to our sum over edge identifications. The 
cycles of the unions of the pairing with each of the first two graphs can be 
interpreted as the vertices in our construction. The two graphs contributing 
to the faces keeps track of the orientation on what may be nonorientable 
surfaces. However, these constructions cannot easily be extended to the case 
where the transposes appear arbitrarily. 

The closed words in [1] may be interpreted as faces, and the condition 
that the graph of a word be a tree when there is one word and a graph with 
exactly one (possibly degenerate) cycle when there are more words can be 
translated into a condition on the number of vertices and hence on the Euler 
characteristic of the resulting surface. 

Approaches such as [5,6,15] may be interpreted in terms of the algebraic 
construction in Section 0J 

Any compact, connected surface without boundary is either a sphere, 
the connected sum of some number of tori, or the connected sum of some 
number of projective planes (for a proof of this classification theorem see, 
for example, [12]). The sphere has Euler characteristic x = 2, the connected 
sum of n tori has x = 2 — 2n, and the connected sum of n projective planes 
has x = 2 — n. Our surface will have the highest Euler characteristic, and 
hence the highest order in N, if it consists of as many spheres as possible. 
The leading terms therefore consist of a disjoint sphere for each face, with 
Euler characteristic 2r, giving it order TV = 1. These terms correspond to 
disjoint noncrossing pairings on each face. 

We will see (Proposition 5.5 ) that the cumulants are the sums of the 
terms corresponding to connected surfaces. (The connected components 
of any surface partition the set of faces, inducing a partition of [r]. For 
any such partition, the contribution to the moment is the product of the 
contribution of each of its components, so the contributions of connected 
components satisfy the moment-cumulant formula Q.) The highest order 
terms contributing to the cumulant k r (Y\, . . . , Y r ) correspond to spheres, of 
order N 2 ~ 2r . 

In the large N limit, only the first cumulant survives. We can also con- 
sider the asymptotic distribution of the Yf~ around their means by examining 
their higher cumulants. The second cumulant is of higher order in TV than 
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the higher cumulants, so asymptotically the distribution is Gaussian (see, for 
example, [13]). The second cumulants can be interpreted as the asymptotic 
covariances of the — E (Yfe). 

4 An Algebraic Expression for the Moments 

Following [15], we now use the methods of [7] to algebraically calculate face 
information from vertex information and edge information (or, in our case 
the dual problem of calculating vertex information from face information and 
edge information). Each end of each edge is labelled. We then construct 
a permutation 7 of these labels which permutes the labels of the edge-ends 
incident to each vertex in a counter-clockwise direction, and a permutation 
7r which swaps the labels of the edge-ends of each edge. The cycles of the 
permutation f 07 then enumerate an edge-end for each edge appearing 
around each face. In this case, however, our surfaces may be nonorientable, 
so we will consider an orientable two-sheeted covering space of each of our 
surfaces. 

An orientable two-sheeted covering space may be defined for any surface 
(see, for example, [11], pages 234-235). This surface is the one experienced 
by someone on the original surface, rather than within it. This construction 
can be performed more simply in our situation, where we have constructed 
our surface as a cell complex. For each face, we construct two preimage 
faces, corresponding to the two orientations, which we will think of as the 
front and the back. For each untwisted edge identification in our original 
surface, we identify the preimage edges of the front faces in the same manner, 
and likewise the preimage edges of the back faces. For each twisted edge 
identification, we identify the preimage edge of each front face with the 
preimage edge of each back face. 

To accommodate this two-sheeted covering, we double the set [2n] on 
which our permutations act. We will use negative integers to represent 
the "back" preimage points, so our set will now be [±2n]. We construct 
a pairing fc which connects the front and back appropriately. We will do 
this by connecting k with I and —k with — / for every pair {k,l} £ ir, then 
"twisting" each pairing by pairing k with —I and —k with I. We will then 
twist each end which corresponds to an X T term. This will result in a total 
of 1 twist if both terms are X terms, 2 twists (equivalent to no twists) if an 
X term is paired with an X T term, and 3 twists (equivalent to 1 twist) if 
two X T terms are paired. 
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Notation 4.1. Let 5 be the permutation on [±2n] which reverses the sign on 
its argument. 

Let (5' be the permutation that takes k to — k and — fc to k if the fcth 
appearance of X in our original expression is transposed and does nothing 
otherwise. 

Conjugating tt by 5 constructs a copy of the pairing on the negative 
integers, and the pairing 5tt5tt consists of both copies of the pairing. To 
put a twist in each pairing, we multiply by 5 again, giving us 55tt5tt = tt5tt. 
We now wish to twist each end corresponding to an X T . We do so by 
conjugating by 5' 

Notation 4.2. Let tt := 5'tt5ttS'. 

We construct permutations taking the place of 7, one acting on the fronts 
(positive integers) and the other acting on the backs (negative integers). 

Notation 4.3. Let 7+ := (1, . . . , 2n\) • • • (2 (n\ + • • • + n r _i) + 1, . . . , 2n) and 
7- := (-1, • • • , -2m) ■ ■ ■ (2 (m + ■ ■ ■ + n r _i) - 1, ... , -2n). 

We would expect the Kreweras complement to be 71-7+ 7Z 1 ; however, 
because of the conventions we are using for numbering the edges and corners, 
we find that instead of negative integer —k, we get 7_ (— k) (while we find the 
correct indices when the integers are positive). To correct this, we conjugate 
by jZ 1 , leaving us with the permutation jZ 1 TT r y+. 

The cycles of this permutation enumerate the subscripts of the matrices 
appearing around each vertex, negative if the matrix must be transposed. 
However, since we are using a two-sheeted covering space of the surface, 
each vertex appears twice. Our application will require us to discard the 
redundant cycles. 

Example 4.4. Let tt = (1, 9) (2, 7) (3, 4) (5, 10) (6, 8). Then if 5' = (1,-1) 
(3, -3) (5, -5) (7, -7) (9, -9), we calculate that tt = (1, -9) (-1, 9) (2, 7) 
(-2, -7) (3, 4) (-3, -4) (5, 10) (-5, -10) (6, -8) (-6, 8), as expected. If we 
have 7+ = (1, 2, 3, 4, 5, 6) (7, 8, 9, 10) and 7 _ = (-1,-2,-3,-4,-5,-6) 
(-7,-8,-9,-10), we calculate that ~{Z l ^ + = (1, 7, -5, -9) (-1, 9, 5, -7) 
(2, 4, 10) (-2, -10, -4) (3) (-3) (6, -8) (-6, 8). 

5 Proof of the Main Formula 

We now present the proof of the main formula. We will need a lemma to 
show the cycles of 7Z 1 7i"7+ do indeed come in front-back pairs as expected, 
so one of each pair may be discarded. In the course of this proof we will show 
some properties of tt which will be useful later. In the more general matrix 
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model with arbitrary transposes (Section |7.3[ ), we will use a permutation 5' 
which might reverse the signs on any k and —k, so the following proofs will 
not depend on these k being odd. 



Lemma 5.1. With permutations as given in Notations 3.1, 4-1 4-% an d 



4-3, if the cycle (ki, . . . , k m ) appears in the cycle decomposition of 7l 1 7T7_|_, 
then the cycle (—k m , . . . , —k\) also appears as a distinct cycle. 

Furthermore, the permutation fc is a pairing with \rr (k)\ = 7r(|A;|). 

Proof. We first show that the cycle (— k m , . . . , — ki) appears. We do so by 
showing that 7Z 1 ^"7+ = 8 (■yZ 1 fcj+) 6; that is, if we reverse all the cycles of 
r yZ 1 fcy+ and change the signs on each integer, we get the same permutation 
back. 

We can see that 5 ("yZ 1 fcj+) 5 = 5yZ[ SStc^SSjS. Since conjugation 
by 5 changes the signs on each integer in the cycle notation of a permutation, 
we can see that Sj^S = jZ and 5j-5 = 7+. 

We can see that fc is self-inverse: fc = (5'tc5tc5') 1 = 5'ic5ic5' , since 5', ic 
and 5 are all involutions. 

We now show that fc is unchanged by conjugation by 5. ir5ir is 5 conju- 
gated by 7r, so its cycle decomposition is that of 5 with each positive integer k 
replaced with I := ic {k). Its cycles are then of the form (I, —k), 1 < k < In. 

Conjugating irbic by 6' reverses the sign on some integers, so the sign on 
either or both integers in the cycle (I, —k) may be reversed. In that case, 
the signs on their negatives in the cycle (—l,k) (which must also appear) 
will also be reversed. So fr = 5'it5tt5' consists of cycles of the form (Z,±fc) 
and {—I, =Ffc), 1, < I < 2n (each cycle will be represented twice here). Since 
k 7^ / for all k, each of these cycles contains two distinct elements, so fc is a 
pairing. We have proven the stated properties of n. 

Conjugating fc by 8 reverses the signs on all integers in the cycle decom- 
position of fc, which changes cycles of the form (I, ±fc) into {—I, ^k) and vice 
versa. So we see that fc is unchanged by conjugation by 5. 

We thus see that 5 (7Z 1 vi"7+) S = 7I" 1 ^7+- We now show that these 
cycles must be distinct. If not, for each i ki and — ki must occur in the 
same cycle. Thus there must be an i and j with i < j such that ki = —kj. 
Then ki + t = —kj—t (by the above arguments) for all integers t. We can 
thus choose our i such that either ki = —ki or ki = —ki + \. Clearly the 
first case cannot occur. In the second case, 7Z 1 ti"7+ (ki) = ki+i = —ki, so 
fc~f + (ki) = 7_ (—ki) = <57+<5 (5 (ki)) = —7+ (ki). However, as we saw above, 
the cycles of fc are of the form (±k, ±7r (k)) (with all four choices of + and 
— possible, depending on k). Since tt is a pairing, it never maps an integer 
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to itself, so 7r cannot map an an integer to another with the same absolute 
value. So it is not possible for n (j + (ki)) to be equal to — r y + (ki). Thus 
ki and — k{ must appear in different cycles, so the cycles (ki, . . . , k m ) and 
(— k m , . . . , — ki) must be distinct. □ 

We can discard the redundant cycles by any method, such as by taking 
only cycles whose integer of lowest absolute value is positive (called partic- 
ular in [9]). 

Notation 5.2. We denote by the set of particular cycles. We will 

also use this symbol to denote the set of integers appearing in these cycles, 
and we will think of the permutation as acting on this set. 

We now prove the main theorem: 
Theorem 5.3. 



a r (Y 1 ,...,Y r )=E(Y 1 ---Y r ] 



N~ 



E ^ z ^ +/2 (Di,...,D 2n ) (5) 

7reP 2 (2n) 



where the notation is as defined above in Notations 3.1, 4-1 4-% 4-3 and 

EH 



Proof. As shown in Lemma 3.2 



a T (Y u ...,Y r ) = N- 



■ E E 

■n-GV 2 (2n) i:[2n]->[7V]:i=io7r 
j:[2n]^[M]:j=joTT 



d (1) d (2) 

'H 7 (l) J2J 7 (2) 



d 



(271-1) 



d? n) 



We will continue to denote -D_& := DT. Likewise, define d) ■ to be the 



«2n-l« 7 (2n-l) J2nJ 7 (2n) 
(-fc) 



,(-*) 



ijth entry of (that is, the jith entry of D^). Then d\- ., ;/ . 

We note that the indices appearing on matrix are ik and l^ha, where 

i represents i or j, depending on which indices appear on that particular 

(k) 

matrix. So the term from which appears in our expression is dl k l , k y If 
we ignore the sign on the subscripts of indices (that is, we define := i^), 
this is equal to d bi ( _ fc)t _ fc . Since 7+ has no effect on negative integers 
and 7_ has no effect on positive integers, both of these expressions are 



equal to d\ 



(*) 



d\ 



(-*) 



S_(fc)S+(fe) - -S_(-*)fcr(-fc)- We can re P lace each term dit\ (k) by 
d[^ k j ±k) t, J+i±k) , where ±k is whichever of k and —k appears in 7Z 1 7T7 + /2, 
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and since exactly one of k and —k appears for each 1 < k < 2n, we can 
choose the superscripts to be the elements of 7Z 1 7i"7+/2. 

We have i = i o tt and j = j o tt, and since \tt (k)\ = tt (\k\) (Lemma 5.1 ), 
Lk = i>a-(k) f° r au L an d &i we can replace the conditions i = i o tt and j = j on 
with i = i o tt and j = j o jr. 

We can now write our expression: 



a r( y 1 ,...,y r ) = iv— E E II <Wc fc r 

TrgT^n) i:[2n}^[N]:i=ion k&~ 1 Z 1 K~ l+ /2 
j:[2n]^[M]:j=joir 

We now show that each of these constraints = l^^) equates a first 
index with a second index in this expression. A second index is of the form 
t 1+ (k) f° r some k G 7Z 1 7i"7 + /2. The index t 1+ (k) is equated with i^ 1+ qa. We 
can see that 7T7+ (k) = 7_ (7Z 1 vi"7+ (&))• The integer 7Z 1 vi"7+ (fe) is the next 
element after k in the cycle of 7Z 1 vi"7+, so it appears in the same cycle and 
hence must also be an element of 7Z 1 7i"7+/2. Then 7_ (7Z 7T7+ (&)) is the 
subscript on the first index on the entry from D -i- r^y Thus the second 
indices are constrained by the values of the first indices, which are otherwise 
unconstrained {tt being a pairing). 

Thus, instead of the constraints i = i o tt and j = j o tt, we can rewrite 
all the second indices as tt^ + (k). This leaves us with terms of the form 

(fc) ^7 (fc) ■ We renumber our indices by letting i' k = i 1 _[k) (legitimate 
since 7_ is a permutation). Our expression is now: 

a r (Y 1} ...,Y r ) = N-"-r E E II d fj 

7reP 2 (2n) i:[2n]-*[N] k£~ / Z 1 *j+/2 
j:[2n]^[M] 

We recognize the term for each tt as the e xpre ssion for the trace along the 
permutation 7r 1 7i"7+/2 as given in Lemma 



2.7 



□ 



We now demonstrate that the cumulants k r (Yi, . . . , Y r ) are given by 
sums over the connected figures, or, in terms of the permutations, those 
where the subgroup of S2 n generated by 7 and it acts transitively on [2n]. 
We demonstrate this by showing that the quantities defined in this manner 
satisfy the moment-cumulant formula 

We note that the elements of a cycle of 7 must be contained in the same 
orbit of (7, tt). 
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Definition 5.4. We say that pairing ir connects cycles ((2 (m + • • • + nt-i) 
+1, . . . , 2 (m + • • • + n k )) and (2 (m + • • • + + 1, . . . , 2 (m + ■ ■ ■ + nj)) 
(which we will refer to as the ith and jth cycle, respectively) if their ele- 
ments appear in the same orbit of (7, 7r). If (7, 7r) acts transitively, we say 
that ir connects the cycles of 7. 

We can construct a partition of [r] where i and j appear in the same 
block when the zth and jth cycles of 7 are contained in the same orbit of 
(7, 7r). We say that ir induces this partition on [r]. 

Proposition 5.5. The cumulants of the Y& are given by 

k r (Y l ,..., Y r ) = J2 Tl lZ^ + /2 i D U • • • > D 2n) 

(7, 7r) acts transitively 

Proof. Define r-linear functions k r by 

~k r (Yi, . . . , Y r ) := ^7 Z ^7 +/2 Pi, • • • , D 2n ) 

(7, 7r) acts transitively 

and for a partition p EV (r) , let 

(^lj • • • j ^r) := J~J &s (iti) • • ■ j ^i s ) ■ 

V={ii,...,i s }ep 

For each p G V(r), each term fc s (Y^, . . . ,Yi a ) corresponding to block 
V = {h, ■ ■ ■ , i s } G p is a sum over all pairings on [2 (n^ + • • ■ + TH a )] con- 
necting the cycles of permutation (1, . . . , 2njj) • • • (2 (n^ + • • • + nj s _ 1 ) , . . . , 

2 (rtij + • • • + n>i s )). If we identify the elements of [2 (n^ -| \- rii s )] with 

the elements of the iith, . . . , z s th cycles of 7 (in order), we can see that this 
permutation is 7 restricted to the elements of these cycles. 

The term k p (Yy, . . . , Y r ) can thus be interpreted as a product of sums 
which, when expanded out, is the sum of terms, each corresponding to a 
choice for each block V G p of a pairing connecting the cycles indexed by 
V, that is, a pairing of [2n] inducing the partition p. Any pairing of [2n] 
inducing the partition p consists of such a choice. Thus k p (Yi, . . . , Y s ) is the 
sum of all terms Tc -1^ , 2 (Di, ■ ■ ■ , D2n) for ir inducing p. 

Since any pairing it must induce some partition p G V (r) , the moment 
a r (Yi, . . . , Y s ) is the sum of the k p (Yl, . . ., Y r ) for all p G P (r). Thus the 
fc s satisfy the moment cumulant formula (|2j). □ 
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6 Asymptotics 



For a £ S n , let # (tt) denote the number of cycles of a (that is, the number 
of orbits in its action on [2n], and for a, r £ S n , let # ((cr, r)) denote the 
number of orbits of the subgroup (a, r). In place of the Euler characteristic, 
we will use the inequality 

# (*) + # (ff _1 r) + # (r) < n + 2# ((a, r)) 

(see [8,15]). 

Since 7Z 1 vr7 + is conjugate to 7_7r 1 vi"7 + 7r 1 = tt^+'JZ , these two per- 
mutations have the same number of cycles, so we find that 

# (*) + # (71^7+) + # (7+7I 1 ) < 4n + 2# ((7+7Z 1 , vr)) . 

We showed in Lemma |5.1| that tt is a pairing, so # (7T) = 2n. We can cal- 
culate that 7+7Z 1 = (1 . . . , 2m) • • • (2 (ni + • • • + n r _i) + 1, . . . , 2n) (— 2ni, 
. . . , —1) • • • (— 2n, . . . , —2 (m + . . . + n r —i)), which has 2r cycles. The sub- 
group (7+7Z ,7r) cannot have more orbits than the subgroup generated by 
7+7Z 1 alone, so it must have at most 2r orbits. The permutation r yZ 1 TT^+/2 
has half as many cycles as 7Z 1 tt7+. So we find 

# (7l 1 7T7+/2) <n + r. 

If we express the moment a r (Yi, . . . , Y r ) in terms of the normalized trace, 
we find that the leading order terms (those for which the equality holds) 
have order iV°. 

When calculating the cumulants, we only consider pairings ir such that 
(7, tt) acts transitively. Topologically, we can see that for such a tt, (j+jZ 1 , tt) 
has at most 2 cycles. Algebraically, for such a tt, we can find a sequence of 
cycles (k,l) of tt connecting the sth and tth cycles of 7 for any s, t E [r]. As 
shown in Lemma 5.1 if cycle (k, I) appears in tt, cycles (k, ±1) and {—k, 
appear in tt, so we can find a sequence of cycles of tt connecting the sth cycle 
of 7+ to the tth cycle of either 7+ or 7_ (or possibly both), and likewise 
with the sth cycle of 7_. Thus the subgroup (7+7Z 1 ,tt) can then have at 
most 2 orbits. Thus 

# (7l 1 7T7 + /2) < n - r + 2. 
The order of k r (Y±, . . . ,Y S ) is then N 2 ~ 2r and 



lim iv 2r - 2 £v(yi,...,y r 



AT- 

V 

N 



£ TviLVW* (D 1 ,...,D 2n ). (6) 



(7, 7r) acts transitively 
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Theorem 6.1. As N — ► oo, the random variable Yk converges weakly to the 
constant random variable with value limjv^oo k\ (Yk). 

Furthermore, the centred and rescaled random variables N {Yk ~ IE(Yfe)) 
converge weakly to a centred multivariate Gaussian random variable with 
covariances limjv^oo N 2 k2 (Yk, Y{). 

Proof. We will use the method of moments in both cases. In each, con- 
vergence in moments implies weak convergence (see, e.g. [17] for a proof of 
the method of moments and, e.g. [2], page 389, for a proof that a Gaussian 
distribution is determined by its moments). 

We can see that the first cumulant of Yk converges to the value given by 
^ as N — > oo, while all higher cumulants converge to zero since they are of 
higher order in jr. These are the cumulants of a constant random variable. 

Centring the random variable Yk does not change the value of any cu- 
mulant k r for r > 1 (since the cumulants are multilinear and any set of 
constant random variables are independent and hence have zero higher cu- 
mulants). After mutliplying by N, the first cumulant is zero, the second is 
limjv^oo N 2 k2 (Yk,Yi) as given in (|6b , and all higher cumulants converge to 
zero as N — ► oo. These are the cumulants of a multivariate Gaussian. □ 

Remark 6.2. As in [15], highest order terms may be described in terms of 
planar diagrams, although we will not do so rigorously. If each face is drawn 
in the plane in one of its orientations, then a pairing corresponding to a 
highest order term is one which may be drawn without crossings while the 
twistedness of each edge identification respects the relative orientations of 
the two faces. 

The asymptotic expected value of a Yk is the same as in the complex 
case, since any twisted edge identification will result in a non-orientable 
surface, which will not be a sphere. The kth moment will have 2 r_1 times 
as many terms as in the complex case, since there are 2 r possible relative 
orientations of the faces, but the bijection is not natural, so unless the 
constant matrices are chosen in a particularly symmetric way, the moment 
will not be 2 r ~ 1 times the value in the complex case. 

7 More General Matrix Models 

7.1 Several Wishart Matrices 

Instead of a single random matrix X, we can consider a family of random 
matrices indexed by a Hilbert space TL: {X g : g £ TC}. We let the ijth entry 
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of Xg be ^hjfij j where the are JV(0, 1) random variables which are 

independent for different choices of i and j and such that E (jlf 1 1 fjjj 
'!■ h) for g,heH. 

Applying the Wick formula creates the same constraints on nonzero 
terms for a given pairing ir. However, the value of the expected value ex- 
pression Q, when it is nonzero, will now be 

n e (/S4??)= n 

{k,l}£ir {k,l}£n 

The term in ([3]) associated to the pairing ir will be weighted by this value, 
so in place of^5|, we get 

a r (Y l ,...,Y r )=N- n - r £ H (g k ,gi)Tr j -i^ +/2 (D 1 ,...,D 2n ). 

neV2(2n) {k,l}eir 

In particular, we can consider several independent Wishart matrices: if 
the pairing connects terms from different matrices, the resulting term will 
be zero. We can colour the edges according to the different matrices, and 
consider only pairings which match colours. 

7.2 q- Wishart Matrices 

We can also consider matrices whose entries are g-commutative Gaussian 
random variables, that is, noncommutative random variables satisfying the 
g-Wick formula, as discussed in [3]. 

Definition 7.1. Let ir £ Vi in) be a pairing. We call a crossing a pair of 
blocks {i, j} , {k, 1} G ir such that i < k < j < I. We denote by cr (tt) the 
number of crossings in ir. 

Definition 7.2. Noncommutative random variables {X g } for some Hil- 
bert space 7i are said to be real q- Gaussian if they satisfy the g-Wick for- 
mula: 

E(X gi ..-X gn )= ? CrW II <9k,9l) 

Again, the expected value expression Q can be replaced by this one, and 
the expression corresponding to pairing ir becomes a weight on the term in 
Q corresponding to ir, so in place of we get 

a r (Yi, ...,Y r ) 

= N- n ~ r 1 CT{n) II (9k,9i)Tr^ +/2 (D 1 ,...,D 2n ). 

7reP 2 (2n) {k,l}eir 
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Remark 7.3. In an expression with only one trace, crossings correspond to 
the geometric notion of crossings. However, in an expression with several 
traces, the number of crossings depends on the order of the terms inside 
each trace, which is not expressed geometrically. 



7.3 Arbitrary Transposes 

We can also generalize our formula to expressions in which X and X T terms 
may appear abitrarily, rather than alternatingly. We let e : [n] — > {— 1, 1} be 
a function encoding the positions of the transposes. We denote by X^ the 
matrix X itself and by X* -1 ) the transpose of X. In place of formula ([I]), we 
can calculate the moments and cumulants of random variables of the form 

Y k = tr (x^ + - n ^+ 1 »D ni+ ... +nk _ 1 ■ ■ ■ X^ + - n ^D ni+ ... +nk 

where the D k matrices are of the appropriate size for the above matrix 

multiplication to be defined. We again let n = n\ H Y n r . It is no longer 

necessary that any of the n/% be even, although if n is odd there will be no 
pairings on [n], so the moments will all be equal to zero. 

Instead of the indices i k and j'^, we sum over all indices 1 < . . . , < 
M and 1 < . . . , < N such that it is the L^L^th. entry of X that is 

taken from its kth. occurrence. We then find that it is the 'o^jl th 
entry of D k which appears, that is, 

d {±k) 

S_(±fe) S+(±fe) 

We now define 8' (k) := e(\k\)k, so it reverses the sign on integer k 
when the \k\th. occurrence of X appears with a transpose, and we construct 
7r using this 5'. Permutations 7, 7_ and 7+ are defined as above. 

The rest of the argument proceeds as above, and we find that 



a r (Y 1 ,...,Y r ) = N-^- r Y, Tv y -^ +/2 (D 1 ,...,D n ). 
Remark 7.4. Since a real Wigner matrix can be written ^ (X + X T ^j, where 



X is a square matrix as defined in Definition 2.2 we can expand to find that 



the nth moment of a compound Wigner matrix can be expressed: 

i N ~^~ r E E ^-,^ +/2 ( Dl ,...,D n ). 

e:[n]-»{-l,l}7reP 2 (n) 

Geometrically, this corresponds to a sum over all possible surface gluings 
with both directions of edge identification possible. 
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Example 7.5. If we wish to calculate the quantity E (TV (XD1XD2)), we use 
5' = e. There is only one pairing on [2]: tt = (1, 2), and tt = (1, —2) (2, — 1). 
We calculate that r yZ 1 TT^ + = (1,-2) (2,-1). Only the cycle (1,-2) is par- 
ticular, so 7Z 1 7i"7+/2 = (1, —2), a permutation on {1, —2}. Thus 

E (Tr (XD 1 XD 2 )) = N~ 2 Tr (D 1 Dj) = 7V _1 tr (Di-Dj) . 
If we let Z := \ {X + X T ) be a real Wigner matrix, then we find that 

E (Tr {ZD 1 ZD 2 )) = ^E (Tr {XD 1 XD 2 )) + ^E (Tr (XD 1 X T D 2 )) + 

^E (Tr (X T D 1 XD 2 )) + ^E (Tr (X T D 1 X T D 2 )) . 
Calculating the last three terms by the same methods, we find that 

E (Tr (ZD 1 ZD 2 )) = ^" 2 Tr (D 1 ) Tr (D 2 ) + ^iV- 2 Tr {D 1 dJ) 
= l -ti{D 1 )ti{D 2 )+ l -N- l ti{D 1 Dl). 

Two sphere terms contribute to the first term on the right hand side, and 
two projective plane terms contribute to the second term. 
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